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Abstract 

The decays 77,7/ — > 7r + 7r~7 are investigated within an approach that combines one- 
loop chiral perturbation theory with a coupled channel Bethe-Salpeter equation which 
satisfies unitarity constraints and generates vector mesons dynamically from composite 
states of two pseudoscalar mesons. It is furthermore shown that the inclusion of the rj as 
a dynamical degree of freedom does not renormalize the Wess-Zumino-Witten term. 
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1 Introduction 



The decay 77 — > 7r + 7r _ 7 is determined entirely by the chiral anomaly, if the quark masses of the 
up, down and strange quarks and the involved four-momenta are sent to zero. The kinematical 
region of the decay, on the other hand, is constrained to 4m^. < (p + +p~) 2 < rrii with p + and p~ 
the four-momenta of the tt + and 7r~, respectively, and thus far from the zero momentum limit 
which is described by the anomalous Wess-Zumino-Witten (WZW) term in the chiral effective 
Lagrangian [1, 2]. 

In the framework of chiral effective field theory higher order contact interactions as well as 
unitarity corrections which arise from loop graphs will be of importance and must be included in 
order to enable a proper description of the experimental decay width and the photon spectrum, 
see, e.g., [3] for a one-loop calculation and [4, 5] for experimental results. 

For the decay t/ — > 7r + 7r~7 contributions from vector meson exchange will dominate the 
amplitude and unitarity effects should be implemented via final state interactions. This is 
clearly beyond the perturbative framework of chiral perturbation theory (ChPT) and requires 
utilization of non-perturbative tools which match onto the results from ChPT. One obvious 
approach would be to employ the vector dominance picture with energy-dependent widths in 
the vector meson propagators. However, this procedure can be shown to be in contradiction to 
the one-loop result of ChPT [6]. 

Another possibility is to postulate an N/ D structure for the decay amplitude which matches 
onto both the one-loop chiral corrections and vector meson dominance in the pertinent limits. 
This approach has been applied successfully in [7], where the decay width and photon spectra 
of both the r\ and rj decay are brought to agreement with experiment [8, 9] and constraints for 
the 7]-i]' mixing angle have been given. 

In the present investigation we will apply an alternative approach which relies solely on 
chiral symmetry and unitarity and has already been employed successfully in the anomalous 
two-photon decays of 7r°, 77 and 7/ [10]. Therein the one-loop contributions from the WZW La- 
grangian were calculated and unitarity corrections beyond one loop were included by employing 
a coupled channel Bethe-Salpeter equation which satisfies unitarity constraints and generates 
vector mesons from composite states of pseudoscalar mesons. Although being similar in spirit 
to the complete vector meson dominance picture, this approach clearly distinguishes between 
the exchange of either one or two vector mesons. It turns out that for the two-photon de- 
cays the exchange of one vector meson is the dominant contribution, whereas the simultaneous 
exchange of two vector mesons is suppressed. This is in contradistinction to complete vector 
meson dominance where the coupling of photons to pseudoscalar mesons is always mediated by 
vector mesons. 

The purpose of the present work is to extend this approach to the decays 77, rj' — > 7r + 7r~7. We 
will first perform a complete one-loop calculation including all counter terms of unnatural parity 
up to chiral order p 6 and with the 7/ as a dynamical degree of freedom, but without employing 
large N c counting rules. As we will show explicitly, the inclusion of the massive 7/ state does not 
lead to a renormalization of the WZW Lagrangian and therefore satisfies constraints from the 
anomalous Ward identities. Unitarity corrections and final state interactions are then appended 
within a coupled channel analysis which can be easily matched to the ChPT result. 

The importance of resonance exchange for these decays can be studied without including 
vector mesons explicitly in the effective Lagrangian. Moreover, we will critically examine the 
issue of 77-77' mixing for these decays. 
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This work is organized as follows. In the next section, the complete one-loop calculation of 
the decays is performed. The inclusion of unitarity corrections beyond one-loop is outlined in 
Sec. 3 and numerical results are presented in Sec. 4 along with a comparison with experimental 
data. Sec. 5 contains our conclusions. 

2 One-loop calculation 

The decays rj, rj' — > 7r + 7r~7 arise from the unnatural parity part of the effective Lagrangian 
which collects the terms that are proportional to the tensor e MI/a 0. Within the effective theory 
the chiral anomalies of the underlying QCD Lagrangian are accounted for by the WZW term 

[1, 2, 11] * 

8 w „(U,v) =J^C„ W = - JLj Jj#) - ^ £ W(U , V) (!) 

where 

W(U,v) = i (U dU ] U dU ] U dU ] v - U ] dU ' U ] dU ' U ] dU v) (2) 

with X = U^dU and for the number of colors we set N c = 3. The matrix valued field U = 
exp{i\/2(j)/ /} contains the Goldstone boson octet (ir,K,r] 8 ) and the singlet field 770, where / is 
the pseudoscalar decay constant in the chiral limit. The expression (...) denotes the trace in 
flavor space and we have displayed only the pieces of the Lagrangian relevant for the present 
work. We utilized, furthermore, the differential form notation 

v = cb%, d = dx^df, (3) 

with the Grassmann variables dx^ which yield the volume element dx^dx u dx a dx 13 = e fJ,ual3 d A x. 
In the second integral Minkowskian space M4 is extended to a five-dimensional manifold M5 
and the U fields are functions on M 5 , cf. [2, 11] for details. The external vector field v = —eQA 
describes the coupling of the photon field A = dx^A^ to the mesons with Q = |diag(2, —1, —1) 
being the charge matrix of the light quarks. 

In addition to the WZW anomaly action there is one more term of unnatural parity at 
fourth chiral order contributing to rj' — > 7r + 7r~7 at leading order and to the 77 decay via mixing 
at next-to-leading order 0{p & ). This gauge invariant term arises due to the extension to the 
U{2>) framework and reads 

d 4 x4? = i W 3 (dU drfdv + dUUU dv) , (4) 

where W 3 is a function of r/ , W 3 (rj / f). The potential W 3 can be expanded in the singlet field 
with coefficients that are not fixed by chiral symmetry, while parity conservation implies 
that W 3 is an odd function of i] . 

At the tree level the contributing diagram is depicted in Fig. 1 with P symbolizing either 
an rj or ?/. Expanding U in Eqs. (1) and (4) in terms of yields the tree level vertex 

3 Note that for our purposes we can safely set the singlet axial vector field (a M ) and the derivative of the 
QCD vacuum angle, 8^9, to zero in S wzw which enables us to work with the renormalization group invariant 
form of the anomaly. 



3 



7T + , p + 




IT ,p 



Figure 1: Tree diagram of the decay P — > 7r + 7r~7, where p + and p~ denote the momenta of the 
outgoing pions. The momentum and the polarization of the photon are indicated by k and e, 
respectively. 



d 4 x£ WZ w = (d<f>d<f>d<f>v) + 

<£x6S = wi l) ^r]o (d(j> d<j)dv) + ... . (5) 



From these terms we can derive the tree level amplitude 

A^\P - vr+vr-7) = -ek&p+p-^-^a™ (6) 

with 

< e£) = ^ a^ = ^j\-l^\ (7) 

and is the coefficient of the leading term in W 3 . 

2.1 One-loop diagrams 

There are four different topologies of one-loop graphs contributing to the 7r + 7r~7 decays. We 
start with the discussion of the tadpole diagram shown in Fig. 2. The pertinent terms from 




7T , P 



Figure 2: Tadpole diagram which contributes to P — > it + it 7. 
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both the WZW and the unnatural parity Lagrangian are 

i\J~2 

d 4 xC WZ w = ~ 24?r 2 ^. 5 ((|</# , dxj>\ [dcf) , d<f>] + #[# , (/# , #]] 

+[[# , , 2 ] + , [d<j> , d(p}d(p}) v) + ... , 

d 4 xC^ = ^ Vo (([<!>, d^, d<j>] - [<j>, [<j>, [d<j>, d<j>}}}) dv) , 

+ 4 3) ^o(##^) + --- • (8) 

The w^-term yields an rj' tadpole which — strictly speaking — spoils the chiral counting scheme, 
as it contributes at fourth chiral order. However, the v( tadpole does not contain any infrared 
physics and can be absorbed completely into the low-energy constants (LECs) of the effective 
Lagrangian, as it is neither a function of the Goldstone boson masses nor of the external 
momenta, it is just a constant. This is consistent with the fact that in infrared regularization 
the rj' tadpole vanishes, whereas the tadpoles for the Goldstone boson octet remain unaltered. 
By defining 

wV=w£»-w<?^/f (9) 

the Vj tadpole is compensated by the second term in Eq. (9) and the chiral counting scheme is 
restored without renormalizing the WZW Lagrangian. The expression A^ symbolizes the finite 
part of the tadpole integral 

d d l i \ ml m m 2 



{2n) d Z 2 - m 2 + is I 16tt 2 ln /i 2 ^ 

™ / finite 

with n being the scale introduced in dimensional regularization. In the present work we are only 
concerned with the finite pieces of the diagrams and neglect the divergent portions throughout. 
Note that there are further contributions from one-loop graphs and counter terms which require 
a renormalization of as will be discussed below. 

The remaining contributions in Eq. (8) to the amplitudes involve only Goldstone boson 
tadpoles and are thus of next-to-leading order 

A (ta d){p ^ 7r+7r - 7) = -ek^p+pjer* _L_ (p™*^ + p ( ; ad) ' K A K ) (11) 



with 



- 5 a(tad),it _ 5 /. / 2 9 „(!)»■ 



Pn - 3^3 ' Pff ~ 3 IV 3 b 

B (tad),K__ 4 B (tad),K _ _5 f [2 lfi 2.,,(l)r 



(12) 



Note that we have replaced the coefficient by the renormalized coupling w^ r which is 
consistent at the one-loop level. 

The loop graph in Fig. 3a has a four-meson vertex from the lowest order Lagrangian £(°+ 2 ) 
of natural parity 

£ (0+2) = y (D^D^U) + . . . = ^_ ( [0 , W , d y }) + . . . , 
D lx U = d l ,U + i[U,v l ,} , (13) 
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7T , p 




7T , P 




(a) (b) 
Figure 3: One-loop diagrams contributing to P — > 7r + 7r~7. 

which yields the amplitudes 

AV(P - 7r + 7r- 7 ) = -efc^p+p-c^^ (tf^AK; + /^'^(m 2 ,; *+_)) (14) 
with s + _ = (p + +p - ) 2 . The integral ii is defined by 

him*;?) = l(^ + (m 2 - ^)G mm ( P 2 ) + ^~ 2 (p 2 ~ 6m 2 )), (15) 

where G is the finite part of the scalar one-loop integral 

d d l i 



{2n) d {I 2 -m 2 + ie)((l - p) 2 -m 2 + ie) 



finite 



16tt 2 



, mm m 2 — m 2 , m 
1 + In — — + - 

/i z p z m 



2 V /\ m jn(p 2 
p 2 



■ artanh 



In — 

m 

\/><mrh{p 2 ) 

m + m) 2 — p 2 



(16) 



A m m(p 2 ) = ((m - mf - p 2 ) ((m + m) 2 - p 2 ) . 
The coefficients /3p^ read 

ofoj.ir 1 ofal.ir /2 2 (l)r 

A? = ^=, /v =y--i67r 2 ^ ; , 

.(a;,* _ 2 _ 1 „ (1)r 



(17) 



The diagram in Fig. 3b, on the other hand, involves the five-meson vertex of the WZW 
action 



— / (S 5 ) = on ^ 2 - (0 ^0 ###) + ... . 



80tt 2 



20tt 2 / £ 



(18) 



7T+, p + 



IT ,p 



7T+, p + 



Figure 4: One-loop diagrams with an rj'ir loop. 

The singlet field rjo does not contribute to the five-meson vertex and the 77' is only involved via 
r)-r)' mixing which is of higher order in this framework. At 0(p 6 ) the loop diagram in Fig. 3b 
thus yields a contribution only to the rj decay 



A 0>)(P vr + vr - 7 ) = -ekw+ppd""** 



1 q(b) A 
-^pPp K 



with 



(b) 



(19) 



(20) 



Two additional one-loop graphs which contribute at sixth chiral order are depicted in Fig. 4 
and an analysis of the 0{p 2 ) four-meson vertex shows that the only contribution to this diagram 
stems from 

£(0+2) = Vl (DprfDHJ) + ... = v^j^l (d^dy) + ... (21) 

where Vi(i] /f) = f 2 /4 + v^i] 2 /f 2 + . . . . The fourth order vertex is given by Eq. (5) and we 
can readily calculate the amplitude 

A (c) {P _> 7r+7r - 7) = _ ek ^ p + p -^_L_ pV^mlm*-, s+7 ) + 7 2 K,mJ; s_ 7 )) (22) 
with s +7 = (p + + k) 2 and s_ 7 = (p~ + k) 2 and 



(c) 



Avf ] I 12 



— 167T W 



2„„( 1 )»" 



(23) 



where we have replaced again by w^ r and neglected two-loop corrections. The finite part 
of the nrj' loop integral I2 is given by 



I 2 (ml,m 2 ,;p 2 ) - (p 2 ~ («V ~ m n ) 2 )(p 2 - (m v + m^) 2 )G m ^, (p 2 ) 

+ (p 2 - m 2 , + ml)K + (p 2 + m% - m\) } 



(24) 
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In order to maintain the matching between the chiral and the loop expansion of the amplitude, 
we have to absorb the contributions of 0(p 4 ) in into the leading order coupling 

2 2 

Eq. (5). Expanding I 2 in terms of and ^f- yields 

v' v' 

(25) 

and the renormalization may be accomplished by modifying Eq. (9) according to 

where we replaced by w<p r in the last term of Eq. (26). After this renormalization, we 
work with the subtracted integral 

I 2 {ml,m 2 v ,;p 2 ) = I 2 {m 2 n , m 2 ,- p 2 ) - J 2 (0, m\,- 0) (27) 

which yields contributions to the amplitude of order p 6 and higher. Thus all one-loop contribu- 
tions of 0(p 4 ) can be absorbed into a redefinition of the low-energy constant w<p and do not 
renormalize the WZW Lagrangian. 

We now turn to the next-to-leading order chiral corrections from the decay constants, rj- 
rf mixing and wave-function renormalization which are readily calculated from the formulae 
given in [12]. Replacing in the tree level result, Eq. (6), the decay constant in the chiral limit 
by the corresponding one-loop expression induces the following 0(p 6 ) corrections to the decay 
amplitudes 

A (f) {P - 7T+7T-7) = -ek^ptp^^-^a^ (^f + 2^f) , (28) 
where the 5F P are defined by F P = f(l + 5F P /f 2 ) with [12] 
6F W = A(3f\2m 2 K + ml) + 4/f ^ "A. - , 
5F V = ^f\2m 2 K + ml) + ^m 2 - ^A K , 

6F„ = 4(2< + ml) {pf> + \p?> - 3/$> + + ffi + 3/$> - /$> - f-^) . (29) 

In the case of the 7/ we employed the QCD renormalization scale invariant expression F v r which 
is derived from the singlet axial-vector matrix element of the 7/ and the (3^ are couplings of 
the natural parity part of the Lagrangian at fourth chiral order, see [10, 12] for details. 
Wave-function renormalization for the external particles and 77-7/ mixing yield 

A< Z >{P - 7r+7r- 7 ) = -ek&p+p^ ^l $> (30) 

with 

^ = ^(^ + 2i2?))+a^< ) . 

= < ee) R?} + ^ (< + 2 ^ 2) ) , (31) 
and the R {2) can be found in [12]. 



W W _ (l)r _ w (3)Xf_ _ V? ( 
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2.2 Counter terms at 0(p 6 ) 

The full set of contributing contact terms of unnatural parity at order p 6 in the U(3) framework 
is presented in App. A. Expanding both the potentials and U in the meson fields one obtains 
(Eqs. (A.3), (A.5)) 



8\/2 

d 4 xC^ = i wf^ — j^- {{x , 0}({^0 d(f) , dv} + 2d(f) dv d<j))) 

+ * ^s 0) —JT~ (x 4>) (# d<p dv) - i r] ({x , , dv}) 

~ i ™io Vo (x) (# d<p dv) 



— iw 



■ (32) 

The first four terms contain the quark mass matrix M. = diag (m, m, m s ), m = m u = ma, which 
enters in the combination x = 2BM. with B — — (0| qq |0) // 2 being the order parameter of the 
spontaneous symmetry violation, whereas the terms %, . . . ,^4 arise from the counter terms 
with five derivatives, since the external vector field counts as order p. The resulting counter 
term contributions to the decay amplitudes read 

A^(P - 7r + 7r- 7 ) = -ek&p+pjeT*-^ (33) 



with 



^ = ^ { " 4«4 0) ™* + 84 0) « - m 2 ) 

+^° 1 ) (mJ + 2s + _ - 2m 2 ) + ™2W " 2m 2 )} , 

^ = 32 7 r 2 V /|{ 8 (-4° ) + **f>)ml + (4tDf + 6^!° ) )(2m 2 , + m 2 ) 

+2wi° 1 ) (mJ + 2s+- - 2m 2 ) + S^m 2 , + «+-) 

+2(^ 2 ) + 3^°3 ) )( S+ _-2m 2 )} . (34) 



The w- '* are understood to be the finite parts of the coupling constants in the Lagrangian of 
sixth chiral order. Within SU (3) ChPT the absorption of divergences from one-loop diagrams 
into the LECs of the 0{p % ) Lagrangian has been discussed in [3]. For the numerical analysis it 
is convenient to decompose the [3p^ as follows, 

/^ = ^ (35) 
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with 



and 



with 



u>(«0 = (-44°) - 2w® - 2wf})ml + &wj?\m 2 K - ml) + w^m 2 q , 
4 s ) = 2^?? + ^° 2 } , (36) 

= 32^(^7^ + + tD«s + _) (37) 



«#> = 2 «>s> + s^s? 

^ = -4(2«4 0) - + wj? + + 3w®)ml + (4^ 0) + Buffi + mj) , 
«#> = 4«>S> + 2<dJ? + Gw«$ + 3w { $ . (38) 

The wty piece in (3^ can be absorbed into 

(1) , -(0) 2 /onA 

™3 = ™3 + 16^V "V ~ * * * ' ^ 39 ^ 

where the ellipsis denotes the renormalization of the one-loop graphs given in Eq. (26) and in the 
following we employ the redefined value for (3^ after renormalization, (3^ = 32tt 2 ^^(w^ 
w^s + -), without changing the notation. 



m) 



2.3 Full one-loop result 

Summing up the contributions up to 0{p & ) we arrive at 

A (l-loop) = A (tree) + j^tad) + jtfa) + ^(b) + ^(c) + ^f/j + + . (4Q) 

The amplitude has the form 

A (l-loo P){ p _ n+n -^ = _ e ^ e ^ p - e ^_J__ (41) 



with 



^P ^P / ^P ^P 



+-^r~ J i( m -; s +-) + -^2— J i( m io s +- 
^p ^p 



+ (/ 2 (mJ, mj; s +7 ) + J 2 (m* , mj,; s_ 7 )) + # } + /^ ct) , (42) 
-^p 

where we have replaced the decay constant in the chiral limit, /, by the decay constants F P 
(P = 77,77') and .F^, Eq. (29), in such a way that the amplitude A. 1 - 1 ' 100 '^ is accompanied by a 
factor 1/(F P F 2 ). Note that in ap ree,) the coefficient has been replaced by w$p r , in order 
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to account for the renormalization of the p A contributions from loop graphs and counter terms. 
For convenience we also show the explicit form of the [3p~ loop ^ 



Pi, 



(l-loop) 



1 1 



- " 2„,,(l)r \ /rm 2 ^2\ V 2 



U/^^-167r^J«-<)-^ 2) 



-3 A^ + h{m\\ s+_) + 2I l (m 2 K ] 



+ 647r 2 (4 m ) +wj s)s +-) f > 



(l-loop) 



167T U> 



2„„( 1 )'- 



P2 



4(2< + ^)(^ + 3^-^3 ) - V / |^ ) ) 



3A, - ^A„ + J^m 2 ; s+ _) + -h{m\- s + _) 



4uJ 2 ) (/ 2 (m 2 , mj, ; s+7 ) + / 2 (m^, mjj,; s_ 7 )) 



4 /2 2 
3 V ir * _m 



2s f A @5,18 2 /2/_(m) . -(a) \ 

w) I fT" ~ ^ J V 3 ^ ^' + ^ ' 



(43) 



where we have used the abbreviations 



P5,18 — P5 + 



2 v "-3 , 
1/^18 



(44) 



This generalizes the result from [3] to the U (3) framework with the rj as an explicit degree of 
freedom. 



3 Generation of resonances 

Resonances, in particular the p(770), play an important role in the decays of r\ and rj' into 
7r + 7r~7. We include them in the same way as for the two-photon decays in [10], i.e. by employ- 
ing a coupled channel Bethe-Salpeter equation (BSE) which satisfies unitarity constraints and 
generates resonances dynamically by an infinite string of meson-meson rescattering processes. 
The scattering potential A is calculated from the contact interactions of the effective Lagrangian 
up to fourth chiral order and iterated in the coupled channel Bethe-Salpeter equation 

T= [t + A-G]- 1 A (45) 

with a modified scalar loop integral 

M , (46) 

where G has been defined in Eq. (16) and a is a real constant chosen in such a way that the 
/j, dependent pieces of G and a compensate. For a diagrammatic illustration of the BSE, see 
Fig. 5. The resulting T matrix describes accurately the experimental phase shifts in both the 
s and p wave channels [10, 13]. 



11 



qurrii 



IT 



q f ,m f 



q f ,m f 





Figure 5: Diagrammatic Bethe-Salpeter equation for meson-meson rescattering. 

In order to implement the non-perturbative summation of loop graphs covered by the BSE 
in the decay processes under consideration, we treat the BSE T matrix as an effective vertex 
for meson-meson scattering. The pertinent diagrams are shown in Fig. 6 and an analysis of the 
remaining loop integrations confirms that — as in the case of the two-photon decays — only the 
p wave part of T contributes. 

There is, however, one difference to the two-photon decays. Whereas for the process P — > 
ry(*)ry(*) the contributions of the coupled channels were at least of two-loop order, the diagrams 
(a) and (c) in Fig. 6 incorporate also one-loop contributions, since T = A at lowest order. In 
order to avoid double counting, we omit those one-loop contributions from the last section which 
are already covered by the coupled channel calculation. Since the coupled channel diagrams 
(a) and (c) exactly reproduce the one-loop results at lowest order, the next-to-leading order 
calculation is still complete. 

The amplitude of Fig. 6a is given by 

A (cc a){p _ 7 (.) 7W) = -eV^ £ -^^Tr a h(ml; T^W) (47) 
with 



"v 



(CCa),K± 
<V 



(m 



K 



,(2) 



(CCa),K°K° _ 



V3 
2 ' 



(CCa),TT d 



= 7, 
1 

" 6 



7 



(CCa),K°K° 



(CCa),K^ 



= -2 



v'(i 18-'' 
2 



v, 







rn 



K 



2 , | A ^5,18 _ 

„.(2) 



(48) 



The symbol in Eq. (47) denotes summation over the meson pairs n + ix~ , K + K~ and K°K° 
and Tp a ^ is the p wave part of the BSE T matrix for scattering of a meson pair a into a meson 
pair b as defined in [10]. The loop integral 1\ is given by 

h(ml;p 2 )=I 1 (ml;p 2 ) + C (f>P 2 



(49) 
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with Ii defined in Eq. (15) and in order to be in better agreement with experiment we make use 
of the freedom to take arbitrary values for the analytic pieces of the integrals which corresponds 
to a specific choice of counter term contributions. 
The diagram in Fig. 6b yields the amplitude 

A (cc b){P _ 7 w 7W) = _ eW + p -^_L_ £' 7 ^A* 

x [f^^iO) A n + f^ K± \0) A K ] , (50) 

where the coefficients ^^ Cb ^' a are given by 

„(CCb),*± (CCb),*± n „(CCb),K± (CCb),K°K° V^3 

IT) Itf ' <T) I Tj 2 ' 



i'r'- K± = -^vw = 2 ^_ {ml _ ml) i ^ _ |_ j (51) 

The diagram in Fig. 6c involves the p wave scattering matrix for the two coupled channels \ttt]) 
and |7r?/), so that the loop integral contains a pion and either an r\ or an r/. 4 The amplitude 
from Fig. 6c and the crossed diagram reads 

A (ccc ){p _ 7 w 7W) = _ eW + p -^_L_ 

i=[/ 2 K,mJ; S+7 ) ff^ + ^ + )( S+7 ) +/ 2 K,mJ; S _ 7 ) if ^")( S _ 7 )] 

16vr 2 ^) [/ 2 («; S+7 ) f^ + ^ + \s +1 ) 

+ 7 2 K,mJ; S _ 7 ) Tf—"'-)( S _ 7 )] J , (52) 

and the integrals J 2 and f 2 have been defined in Eq. (24) and Eq. (27), respectively. 

We also consider the diagram with two coupled channels depicted in Fig. 6d. The pertinent 
amplitude reads 



1 I 1 

x -< 
2 



+ 



A (*cc ){p ^ 7 (.) 7W) = -eh^plp^^^r^ 

^ a,b 

x i 1 {m 2 a ; 8+ -) f^is^) A^f^iO) A^ + f^^O) A*] . (53) 



4 The (ir, rj) loop stems from the 0(p 4 ) part of the scattering potential A and did thus not appear in the 
next-to-leading order calculation, where only the four-meson vertex of second chiral order was considered. 
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with coefficients ^ cc ^ a,b symmetric under a <-> b 

(2CC),tt ± ,K ± _ (2CC),tt ± ,K°K° _ _\ (2CC),K ± ,K°K° _ 

it] ir) 2 ^ 4 ' 

(2CC),tt ± ,K ± _ _ (2CC),tt ± ,K°K° _ 1 (2 CC),^ ,K°K° 
Irj' Tjj' 2 7r)' 

= - "4) - C I (54) 



oV""K ""tt7 I - p 2 (2) 



and zero otherwise. 



4 Numerical results 

The spectra of the decays rj, rj' — > 7r + 7r~7 have been measured with high statistics [4, 5, 9, 8] 
showing that the rj' decay is clearly dominated by the p resonance. Therefore, the next-to- 
leading order calculation of Sec. 2 is insufficient to describe the rj' decay and we will discuss 
the one-loop results only for the r] decay. As a first estimate we use a set of parameters which 
is consistent with previous one-loop calculations in U(3) ChPT [10, 12, 15] 

F n = 92.4 MeV , F v = l.3F n , F v > = l.lF n , 
4 0) = 1.4 • 10- 3 , = 1.2^/4, (55) 

while neglecting all other LECs including the coupling constants of the unnatural parity La- 
grangian of sixth chiral order. For the regularization scale p we use 1 GeV. The decay width 
turns out to be r n = 14.48 eV which is by a factor of four smaller than the experimental value 
quoted by the Particle Data Group [14] 

T,, = (56.1 ± 5.4) eV . (56) 

A mixing parameter of v!p = 0.5F^/4 corresponds to an even smaller decay width of T v = 
8.59 eV, whereas a value of = 1.5F%/4 results in r n = 17.47 eV. Hence, for realistic values 
of the r]-r]' mixing parameter it is not possible to match the experimental decay width 
without including the contact terms of 0(p 6 ). 

However, both the photon spectrum and the partial width of the decay rj — > 7r + 7r~7 can 
be reproduced at next-to-leading order, if 0{p & ) counter terms are included. But the choice of 
parameters is not unique, hence it is not possible to constrain the mixing parameter from 
a fit to the experimental data on the 77 decay, since variations in can be compensated by 
tuning the counter terms of unnatural parity; larger positive values of e.g., require smaller 
counter term contributions. For consistency with previous one-loop calculations in U (3) ChPT 
[10, 12, 15] we prefer to use a mixing parameter of = 1.2F%/4. At the one-loop level the 
constants w^ r and appear in a linear combination and thus cannot be fixed separately. 
In order to work with a minimum set of parameters, we set in the one-loop calculation w^ r 
to zero and choose the combinations of counter terms Eq. (36) to be = —0.31 ■ 10~ 3 and 
4 S) = 12.9 • 10~ 3 GeV -2 . 
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Figure 7: Photon spectrum resulting from the next-to-leading order calculation with p 6 counter 
terms (dashed) and from the full calculation including the coupled channels (solid). The dotted 
line corresponds to the simplest gauge invariant amplitude A oc k^p^p^e^" 13 . The data are 
taken from [4]. 

The one-loop result is compared with experimental data from [4] in Fig. 7 and shows good 
agreement after accounting for the detection efficiencies. Good agreement with the similar 
experiment [5] is also achieved. It should be mentioned that the counter term contribution 

— (s) 

yields an important contribution to the decay amplitude and dominates the one-loop 
corrections. This fact is a reflection of the tail of the p resonance in the perturbative expansion 
[3]- 

We now turn to the discussion of the coupled channel calculation. The computation includes 
the tree level graphs and all next-to-leading order corrections. In order to avoid double count- 
ing, the one-loop contributions which are already covered by the coupled channel diagrams in 
Figs. 6a, c have been subtracted. The results are plotted in Figs. 7, 8 and compared with the 
experimental data. The agreement with the data is very good and shows that the important 
degrees of freedom are incorporated in our model. For consistency with previous calculations in 
this framework [10, 16], we have set the mixing parameter and the renormalized coupling 
constant w<p r from the unnatural parity 0(p A ) Lagrangian to zero. We furthermore neglect 

(2) 

the unknown LEC v\ , in order to work with a minimal set of free parameters. For the fit of 
our results to the central experimental values of the partial decay widths, T v = 56.1 eV and 
T v i = 59.6 keV [14], respectively, we employ for the subtraction constant C n in the outer loop 
integral with pions (Eq. (49)) the value C n (/j, = 1 GeV) = —1.42 • 10~ 2 which is comparable 
in size with the one used in [10] (C n = — l/(67r 2 ) —1.69 ■ 10~ 2 ), while setting those for the 
other loops to zero. Since different sets of counter term contributions are summarized in C n , 
the values of these subtraction constants may differ for the decays into or two photons. 
The counter term contributions from the Lagrangian of sixth chiral order are needed to bring 
our results to better agreement with the shape of the experimental spectra. The values for the 
combinations of LECs in Eqs. (36) and (38) are 

4 m) = -2.21 • lO" 3 , wi) s) = -9-50 • 10~ 3 GeV" 2 , 
= -6.90 • 10- 3 , w^= 1.00 -10- 3 GeV" 2 . 

It should be emphasized that the choice of parameters is not unique, since for rj, rf — > 7r + 7r~7 
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Figure 8: Invariant mass spectrum of the 7r + 7r~ system resulting from the full calculation. The 
curves are normalized to the integral of the experimental distribution. Data: left diagram [8], 
right diagram [9]. 



variations in one of the parameters may be compensated by the other ones. However, with the 
choice in Eq. (57) the remaining set of parameters is in agreement with previous work [10, 16]. 

In general, the values of the counter terms in the non-perturbative coupled channels ap- 
proach differ from those in the one-loop calculation as already observed for the two-photon 
decays [10] and the hadronic decays of 77 and rj' [16], since in the loopwise expansion the ef- 
fects of resonances are hidden in the LECs, whereas they are generated dynamically in the 
non-perturbative approach. 

In Fig. 9 we show the dependence of our results on the r]-r]' mixing parameter In both 
decays the heights of the spectra are reduced for increasing values of yielding smaller decay 
widths. As for the two-photon and hadronic decays, ~ is thus in better agreement with 
the data. 

The results for different values of the LEC w^ r are depicted in Fig. 10. When keeping 
= 0, the i] decay is not affected significantly. Similar to the two-photon decays negative 
values for w<p r increase the p peak while positive values decrease it. Changes in w^ r could 
in principle be compensated by altering the subtraction constant C n in the loop integral I\. 
The influence of v± on our results is rather small. Variations within a natural range for a 
coupling constant in the 0(p 2 ) Lagrangian 5 , —3.0 • 10~ 3 GeV 2 , . . . ,3.0 • 10~ 3 GeV 2 , yield only 
slight corrections. 

Finally, we compare the contributions of the different coupled channel diagrams in Fig. 6. 
From the plots in Fig. 11 we see that the amplitude for the pure 7r + 7r~ final state interac- 
tion, J\( CCa ) (Eq. (47)), furnishes by far the dominant part. The diagrams which involve the 
five-meson vertex from the WZW Lagrangian (Fig. 6b, d) yield small corrections. This is in 
contradistinction to complete Vector Meson Dominance, where only these two diagrams are 
present. The contributions from the coupled channel diagram in Fig. 6c are almost negligible. 

5 The coefficients of (d^Wd^U) and {W X + X^U) arc / 2 /4 w F%/A w 2.22 • lO" 3 GcV 2 , for comparison. 
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Figure 10: Dependence on the LEC w^ r : w^ r = —2.0 x 10~ 3 (dashed), w^ r = (solid), 
*4 1)r = 2.0 x 10~ 3 (dotted). 
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Figure 11: Comparison of the different coupled channel contributions in Fig. 6: Full result 
(solid), result without A (CCb) (dashed), without A (CCc) (dotted), without A (CCd) (dot-dashed). 
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5 Conclusions 



In this work, we have calculated the anomalous decays rj, rf — > within an approach 

that combines ChPT with a non-perturbative scheme based on coupled channels. This method 
satisfies unitarity constraints and generates vector mesons from composed states of pseudoscalar 
mesons without including them explicitly in the effective Lagrangian. It had recently been 
applied in the anomalous sector for the two-photon decays of tt°, 77 and rj' and is now extended 
to the decays into 2^7. 

We first performed a full one-loop calculation within the framework of chiral perturbation 
theory and without imposing large iV c counting rules. The presence of the massive state 77' 
spoils the strict chiral counting scheme such that loop diagrams with an 77' also contribute at 
fourth chiral order, i.e. at the same order as tree level contributions from the WZW Lagrangian. 
We have shown explicitly that these contributions can be absorbed into a non-anomalous con- 
tact interaction of unnatural parity, but they do not renormalize the WZW term so that the 
anomalous Ward identities are preserved. 

If the counter terms of unnatural parity at subleading chiral order p 6 are omitted in the 
one-loop calculation, the results are in disagreement with the experimental decay widths. Tak- 
ing 77-77' mixing into account ameliorates the situation slightly for the 77 decay, but is still in 
contradiction to experiment. As a matter of fact, it is not possible to constrain 77-77' mixing 
from the decay 77 — > 7r + 7r~7. Moreover, one should keep in mind that mixing effects are of sixth 
chiral order, if large N c counting rules are not employed, and enter at the same order as the 
neglected counter terms. The inclusion of contact interactions of sixth chiral order improves 
the situation for the 77 decay, whereas the results for the 77' decay cannot be brought to agree- 
ment with experiment by adjusting the counter terms of sixth chiral order. This is due to the 
presence of vector mesons, which dominate the 77' decay. While for the 77 decay only effects of 
the tail of the resonances contribute which can be treated perturbatively by absorbing them 
into the couplings of the effective Lagrangian, one must include unitarity effects via final state 
interactions in the 77' decay. 

The inclusion of the coupled channel formalism provides a framework which reproduces both 
the decay widths and the experimental dipion invariant-mass spectra of the 77 and 77' decays 
while matching onto the results of ChPT. Our choice of parameters is consistent with that of 
the two-photon decays as discussed in [10] and provides a non-trivial check for our approach. It 
furthermore confirms that the inclusion of unitarity effects via the Bethe-Salpeter equation can 
be accomplished in the anomalous sector of QCD in a similar way as for the hadronic decays 
of 77 and 77' [16]. Finally, we have shown that this method does not renormalize the WZW term 
and hence satisfies constraints from anomalous Ward identities. 
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A 0(p Q ) Contact Terms in P — > tt + tt~^ 

In this section we discuss the 0(p 6 ) counter terms which contribute to 77 and 77' decays into 
7r + 7r - 7. We use the notation of [10] for the covariant derivative D^U and the field strength 
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tensors R^ u and L^ u of the right- and left-handed external fields, respectively, and make use of 
the following abbreviations: 

M =W X + X j U, N =W X -X i U, (A.l) 

C M = WD^U , Ep, = WDpDM - (D^DMYU . 

These are the building blocks for the construction of the counter terms. Essentially, there 
are two types of counter terms; either they include the mass matrix Ai = diag(m, rh, m s ) 
or they contain two additional derivatives resulting in contributions to the decay amplitude 
proportional to s+_ which can be expressed in terms of the photon energy u. In the SU{3) 
framework the entire set of 0(p 6 ) terms of unnatural parity has been presented in [17], but the 
inclusion of the singlet field r] induces additional structures. The mass terms may be written 
as 



+W 8 ((MC,) (c v Q«p) + (AO (P^C a C p )) 

+W 9 [(MiQ^Ca + C a Q, v )) + (iV(P M ,C a - C Q P^))) (C p ) 

+W W (M) (Q, v C a ) (C» } , (A.2) 

where the last two terms arise from the extension to U(3). Expansion in the meson fields yields 
in the differential form notation of [11] 

d 4 xC% = iwf ) ^({ X A}{{d<pd<P,dv} + 2d<pdvd<P)) 



_ (0 ) 32_\/2 (d(pd(pdv) - iwf } ^? 



+ iw 8 —p~ M) (d(f)d(f)dv) - iw\ —^-^o ({x , d<f>}[d(f> , dv]) 
16^6 

T 3 



and we have integrated by parts in order to simplify the structures. The terms with five 
derivatives and one external field are 

Cf = e^{w n (P, u (E\CpC x - C X C P E\)) 

+W 12 (P^{E\C x C fi - C P C X E\)) 

+W 13 (P^(E X a C x - C X E X S) (Cp) 

+W U (P, V {E\C P - C P E\)) (C x ) } (A.4) 
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with Wi 3 an d W14 being only present in the U (3) framework. Expansion in yields the vertices 



d 4 xC ( ^ ct = - i w ( {? < (d x d(f) dcj) <9 A + d X (j) d<p d x d<p) dv) 

- * w>i2 ( ^ d( ^ d ^ d ^ + W d ^ dv ) 



P 

iw^^(d<t>){[d x d<t>,d x <t>\dv) 



- * S> ^ <[9 A # , rf W > . (A.5) 
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